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Abstract 

We propose to investigate the properties of nuclear matter by measuring the elastic 
scattering of J/ip on nuclei with high precision. The J/ip mesons are produced from the 
photons emitted in high energy electron-proton or electron-nucleus scattering in the low-x 
region. The measurement could be performed at the future ENC, ETC or LHeC facilities. 



1 Introduction 

Nuclei have fascinating properties which are very interesting for basic and applied science. A 
better understanding of nuclei could lead to efficient energy production and other industrial 
uses. The basic models of the nuclei - the liquid drop and the shell model - give a simple but 
refined understanding of many nuclear properties. They are built on the fact that protons 
and neutron in nuclei behave like an incompressible fiuid and are arranged in highly ordered 
quantum states, or shells. A familiar but very puzzling property of nuclear matter is the 
growth of its volume proportionally to the number of nucleons. Although often explained 
by the Fermi nature of the nucleons it is mysterious when viewed from the particle physics 
perspective because nucleons consist of quarks and gluons which are point-like and nearly 
massles^. We note that the sizes of atoms are nearly independent of the charge number Z. 

One of the main reasons why, after almost 80 years of investigations of nuclear structure, 
some of the basic properties are so poorly understood is a lack of proper tools to view inside 
nuclei. In the past, the most direct information about nuclear structure was obtained by 
scattering electron beams off nuclei. The electron is a very good probe because it penetrates 
into the nuclear interior without being absorbed. Unfortunately, electrons can only see the 
electric charge distributions. They are not sensitive to the distribution of strong fields which 
keep nuclei together. Another important source of information comes from the scattering of 
low energy protons on nuclei. The protons sense the full matter distribution but they are 
themselves complicated objects whose interactions with matter are not very well known. Our 

^The size of the hadrons made of quarks and gluons is also not understood. 
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understanding of strong interactions is currently only robust at distances much smaller than 
the proton size, where perturbative QCD is applicable. Although Lattice Gauge Theory has 
firmly established that QCD is the correct theory of strong interactions at large distances, 
its applications to hadronic interactions are only now being developed [1-3]. 

In recent years HERA has shown that, in deep inelastic electron-proton scattering at 
small X, up to 20% of events are of diffractive origiii^. This is comparable to the rate of the 
elastic reactions in hadron-hadron scattering at high energy. In hadronic reactions the rate 
of elastic and inelastic events is connected by the optical theorem. Before HERA, the optical 
theorem was not expected to play any role in electron-proton scattering because the elastic 
ep component is miniscule. The abundance of the diffractive (i.e. quasi elastic) processes 
seen at HERA indicates therefore that some intermediate states have to be formed to which 
the optical theorem can be applied. The theoretical framework which describes naturally the 
formation of such intermediate states is the dipole picture. Dipoles are small quark- ant iquark 
pairs which interact by color exchange. In contrast to hadron scattering, where the projec- 
tiles are complicated, the color dipoles are simple objects. Their creation, interaction and 
annihilation has a simple description in perturbative QCD combined with QED. A multi- 
tude of HERA measurements such as inclusive cross-sections, charm production, diffraction, 
jet production, vector-meson production and DVCS cross sections are accurately explained 
within the dipole picture [4-13] . 

A particular final state into which a dipole transforms acts as an analyzer. For example, 
jets select very small dipoles, J/ ip production selects charm anti-charm dipoles, p mesons uu 
or dd dipoles and so on. For nuclear structure investigations the most interesting process is 
the scattering of small dipoles. The smallness of the dipole assures that the interaction with 
the nucleon is well described by perturbative QCD. In the low- a; region the small dipoles 
interact with the nucleon by gluon exchange only. The color exchange is of short range in 
contrast to the interaction of electromagnetic dipoles. Of particular interest is the elastic 
scattering of color dipoles because the transverse deflection of the dipole measures the spatial 
distribution of the gluonic matter. 

In this paper we will concentrate on the measurement of exclusive J/ip vector meson 
photoproduction. The J/ip meson is a bound state of a charm quark antiquark pair and 
therefore the corresponding dipole is naturally small. The cross section is relatively high 
because we take photoproduction. In nuclei, the small charmed dipole scatters presumably 
on the individual nucleons. In spite of the high energies involved, the nucleus will frequently 
remain intact because the large absorption cross section together with the optical theorem 
assures that the scattering process has to be coherent in around 15% of cases [13,21]. Since 
the charmed quark dipole interacts almost completely via two gluon exchange with matter 
the deflection of the J/tp measures directly the intensity and the spatial distribution 
of the strong field which keeps the nucleus together. This will allow to see precisely the 
structure of gluonic fields which was never seen before. 

The paper is organized as follows; in Section 2 we give an overview of the main properties 
of the dipole interaction and summarize in some detail the existing experimental results. In 



^Here diffraction means tliat the proton stays intact. 
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Section 3 we discuss the properties of J/ip as a. probe of proton and nuclei. In Section 4 we 
discuss the J/ ip production cross section within potential detector acceptances and in Section 
5 we discuss the detector requirements necessary to perform a precise pt measurement and we 
discuss the experimental requirements necessary to assure the elastic signature. In Section 6 
we summarize the results and in Section 7 we conclude. In the Appendix we give a derivation 
of the dipole representation in a simple case of elastic photon-proton scattering. 



2 Dipole description of DIS reactions 



Dipole interactions are, in a way, as fundamental as Rutherford scattering. In Rutherford 
scattering the incoming electron emits a virtual photon which scans the charge distribution 
of the nuclear target. In elastic dipole scattering the incoming electron emits a virtual photon 
which turns into a small quark antiquark pair. In leading order QCD, at large virtuality 
scales and high energies, the qq pair interacts elastically with the nucleon by exchanging two 
gluons with large transverse momenta, / and / + A. The transverse momenta of the gluons 
cannot be directly observed, however their difference. A, is measurable in elastic vector 
meson scattering because it is equal to the difference between the transverse momentum of 
the incoming virtual photon and the final vector meson. The two gluon interaction leads 
in QCD (together with QED) to the dipole representation. In the Appendix we derive the 
dipole representation from the Feynman diagrams in a simple case of virtual photon-proton 
elastic scattering, when A = 0. In the following we discuss the main properties of the dipole 
representation. 




Figure 1: Elastic scattering of a. J/ip meson on a proton or ion in the dipole representation. 
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2.1 Dipole representation 



The dipole interaction proceeds in three stages: first the incoming virtual photon fluctuates 
into a quark- antiquark pair, then the qq pair elastically scatters on the target, and finally 
the qq pair recombines to form a final state. The creation, scattering and recombination of 
the dipole occur in different space time regions because the life time r of the qq fluctuation 
is very long. In the proton rest frame, r ~ l/rripX. Even at a; = 10~^ this correspond to the 
distance of ~ 20 fm, which is larger than the size of all nuclei, i? ~ 1.2 ■ A^^^ fm. In addition, 
due to high energies, the transverse positions of the quark and antiquark do not change in 
the scattering process so that the dipole is not changing its size (see Appendix). 

The amplitude for the elastic photon-proton scattering, 7*p — > 7*p, A^*^{x,Q, A), is 
therefore the product of amplitudes of these three subprocesses integrated over the dipole 
variables r and z: 

A'^Pix, Q,A) = J2j d'r ^*{r, z, Q) A,,{x, r. A) M/(r, z, Q), (1) 

where \E'(r, z, Q) denotes the amplitude for the incoming virtual photon, with virtuality Q, 
to fluctuate into a quark-ant iquark dipole with flavour The wave function is determined 
from light cone perturbation theory to leading order in the fermionic charge, see Appendix. 
The amplitude for the qq to recombine to a virtual photon is \E'*(r, z, Q). Aqq{x, r, A) is the 
elementary amplitude for the scattering of a dipole of size r on the proton, A denotes the 
transverse momentum transfered from the dipole to the target, and x is the Bjorken variable. 
The sum should be taken over all quarks flavours /, including charm. For charmed quark 
the deflnition of the Bjorken x = xb should be replaced by a; = a;s(l -|- Am'l^/Q'^) to take 
into account the charm threshold effects. For similar reasons, in the case of vector meson 
production, x = + My/Q"^). 

The elementary elastic amplitude Aqq is deflned such that the elastic differential cross 
section for the qq pair scattering on the proton is 

^ 'A,(x,r,A)|% (2) 



dt IGtt 



where t = — A^. The notation follows the conventions of ref. [9,18]. Elastic dipole scattering 
is connected by the optical theorem to the total cross section for photon-proton or photon- 
nucleus scattering. To evaluate the connections between the total cross section and various 
diffractive reactions it is convenient to work in coordinate space and deflne the S-matrix 
element at a particular impact parameter b 

Aqq{x, r. A) = j d^b e-'^-^ Aqq{x, T, 6) = i j d^b e"'^^'^ 2 [1 - S{x, r, b)] . (3) 

This corresponds to the intuitive notion of impact parameter when the dipole size is small 
compared to the size of the proton. The optical theorem then connects the total cross section 



^We suppress here references to the photon hehcities for simphcity. 



4 



for the qq pair scattering on the proton to the imaginary part of the forward scattering 
amphtude: 

.«(.,.)^ImA,-(^,.A^0)^/d^S2|l-ReS(.,..)]. (4) 

The integration over b of the S-matrix element motivates the definition of the differential 
dipole cross section as 

^ = 2[1-Re^(a:,r,6)]. (5) 

The total cross section for 7*p scattering, or equivalently F2, is obtained, using ([T]) and 
dl]), by integrating the dipole cross section with the photon wave functions: 



a^;iix,Q) = lmAl\l{x,Q,A = 0) = ^ d¥^'g(vI/*vI/)^_^ r), 



(6) 



where t or ^ denotes the transverse or longitudinal polarization of the incoming photon. For 
completeness we also give the relation between the F2 structure function and the Yp cross 
section at small-x. 

Elastic vector-meson production appears in a similarly transparent way. The amplitude 
is given by 

A^*p^pviA) = J cfr J ^ J d^b-^*y-^ exp{-ib ■ A)2[l~ S{b)]. (7) 

Assuming that the S-matrix element is predominantly real, we may substitute 2[1 — S{b)] 
with daqq/d%. Then, the elastic cross section is 



dt 167r 



dh I — ! d26^t^^exp(-z6- A)'^'''^^" 



47r / ^ ^ d% 



The equations ([6]) and ([H]) determine the total DIS cross sections and the exclusive diffractive 
vector- meson production cross section. Fig. [1] shows a diagram of the dipole scattering when 
the final state is a. J/ip meson. 



2.2 Dipole cross sections 

The universal dipole cross section do"qg/d^6 contains all the information about the gluon 
content of the target and the QCD evolution of the gluon density. Its particular form 
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depends on the evolution schema. In the DGLAP formahsm, the dipole cross section for 
small dipoles is given bjfl, 

^ = yr2«,(/.^)x^(x,/i2)r(6). (9) 

Here xg{x,fi^) is the gluon density and yU^ denotes the evolution scale. The starting scale 
is denoted by fi^. The scale is given by the inverse of the dipole size and is usually taken 
as /i^ = + /ig, to stabilize the behavior for large dipoles. The function T{b) is the 
transverse profile of the gluon density of the proton. It is convenient to write the gluon 
density in a semi-factorized form, xg{x, fj,'^)T{b,x). Because the x dependence of T{b) is 
weak it is omitted in the following discussion. At every x the transverse profile is normalized 
to 1, J (P'hT{h) = 1. The parameters of the gluon density are determined from the fit to the 
total inclusive DIS cross section using eq. [61 The fit is shown in Fig. [2l The gluon density at 
a scale /Xq is parametrized by xg{x, ^l) = c ■ x~^. The three free parameters are c, A and fi^. 
The prediction of the model, at low values, also depend on the assumption on the quark 
masses, rrif. Other schemas like the CGC evolution or Regge parametrizations lead also to 
a successful description of the data [19,20,31]. More details can be found in [18,19]. 

The transverse profile is assumed to have a two dimensional Gaussian form: 

Tib) = -^expi-by2BG). (10) 

since the measured ^-distributions were found to be well described by the exponentially 
falling distribution, daylj/dt oc exp(— _B£)|t|), see Fig. [^. The coefficient -B_d is not exactly 
equal to Bq because the size of the vector meson is not negligible (see the discussion in 
Section 2.4). In addition, the observed coefficient Br, shows some x dependence (e.g. in J /ip 
photoproduction) and so Bq should be x dependent (see the discussion in Section 3.1). We 
note also that the fit with the dipole form factor, da/dt oc 1/(1 — t/M'^Y, which is usually 
used to parametrize the proton form factors does not describe data appropriately [25] . 

The dipole cross section, eq. [9], is propotional to the dipole area, r^, and the gluon density 
seen by the dipole, xg^Xjfi"^). It is therefore a measure of the gluon density. When either 
the dipole size or the gluon density is not small it is convenient to use the eikonalized form 
of the dipole cross section to take into account possible saturation effects, 

-2(1- ,,^^_ -'r'<^s{^^')^9{x,^')Tib) \ ^ /i _ eM--)) ■ (H) 



d% \ 2-3 7 V 2 ^ 

Here Q denotes the opacity which is equal to the right side of eq. M The above formula is 
called the Glauber- Mueller dipole cross section. A diffractive cross section of this type was 
used around 50 years ago to study the diffractive dissociation of deuterons by Glauber [26] 
and reintroduced by Mueller [27] to describe dipole scattering in deep inelastic processes. 
In contrast to Glauber scattering, in DIS the functional form of the opacity is known from 
QCD. 



■^Thc derivation of the dipole cross section is given in the Appendix 
The Fourier transform of a gaussian in 6 is a gaussian in A corresponding to an exponential in t. 
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Figure 2: The 7*p cross section as a function of W"^. The sohd hues are showing the dipole 
fit to the data [17] in the low-x region, x < 10~^. 
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2.3 DVCS and VM production 



The elastic cross section, eq. [H was derived under the assumption that the dipole size is much 
smaller than the proton size. For dipoles with a size r the explicit QCD calculation [45] shows 
that A conjugates to b+ (— 2;)r. Therefore, the cross section is not only sensitive to the proton 
impact factor, b, but also to the dipole size, r. The modified dipole cross section is then 



dt 167r 



j j ^ j ^i'&^^^exp(-z[6+(l-2)rl- A)^ 



d% 



(12) 



This cross section can also be used to describe the DVCS process for which the final state 
consists of the scattered electron, proton and a real photon. The wave function for the 
outgoing state is just the amplitude \E'*(Q^ = 0) for the real photons. Equation \n\ then 
gives, after small correction for the real parl|§ and skewednes^ are applied, an absolute 
prediction for the DVCS process. The prediction is shown along with HERA data in Fig |H 
The figure shows an impressive agreement between the dipole model predictions and the 
data if we reahze that two very different processes are described with the same amplitudes; 
the average event which contributes to the total cross section has around 40 particles in the 
final state, while the DVCS event has just a proton, a photon and an electron. The wave 
function of the final state, \E'*(Q^ = 0), is very different than that of the incoming virtual 
photon. And yet all the distributions and the absolute event rates of the DVCS process are 
properly described. This means that the net effects of possible higher order QCD radiative 



^The derivation of the cross section for exclusive vector meson production or DVCS relies on the assump- 
tion that the scattering amplitude is purely imaginary. The real part of the amplitude can be accounted for 
by multiplying the exclusive cross section by a small correction factor given in ref. [11]. 

^For vector meson production or DVCS one should use the off-diagonal gluon densities, since here the 
two gluons carry different fractions x and x' of the proton momentum. This effect can be taken into account 
by multiplying the gluon distribution xg{x,fi^) by a correction factor, Rg given in [29]. 
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Figure 4: The total DVCS cross sections cr as a function of (left) and W (right) [32] 
compared to predictions from the b-Sat and b-CGC models [19]. 



corrections [28] or possible kinematical corrections [30] , have to be small or are consistently 
absorbed into the gluon density. 

Successful comparisons of the inclusive diffractive cross sections with the dipole model 
predictions were performed in several investigations [7,8,13,31]. Inclusive diffraction denotes 
the sum of all processes in which the proton remains intact. Since the summation is over a 
complete set of states the cross section depends only on the photon wave functions. However, 
the dipole cross section in this case has to take into account also the possible qqg Fock state 
which makes the comparison less clean than in the DVCS case. Nevertheless, the comparison 
of the cross sections with the predictions works again very well [13]. The same is also true 
for the comparison with the diffractive jet cross sections [37]. 

Further support for the validity of the dipole picture comes from the very good agreement 
of the predicted vector meson cross sections with the data, shown in Fig. [51 The vector meson 
wave functions are constructed according to a standard procedure developed in [43,44]. 
They are educated guesses based on general considerations described in [11], [18]. The 
only phenomenological input which is necessary is obtained from the measured electronic 
decay width of the vector meson. The very good agreement between data and dipole model 
predictions means that the vector meson wave functions were properly estimated. This is 
further supported by the excellent description of ai/dT ratio as a function of for the J/ip 
and vector mesons [18]. In case of p mesons this ratio is not so well described. This can 
be attributed to the complicated dynamic of the p decay into pions. 
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Figure 5: The total cross section a vs. (Q^ + My) for exclusive Jjil) [25,33], [35] and 
p [36] meson production compared to predictions from the b-Sat and b-CGC models using 
the "boosted Gaussian" vector meson wave function [19]. 

2.4 The ^-distributions 

At HERA, the distributions of the square of the momentum transfer, t = — A^, have been 
measured in exclusive vector meson and DVCS processes. The value of t is usually determined 
from the pt imbalance of the final state particles seen in the detector. The t distributions 
have also been measured using forward proton spectrometers for all diffractive processes. 
The forward spectrometers have however a relatively small acceptance. 

All measured t distributions are very well described by the exponential dependence 
da/dt oc exp(— The t range used in the fit to the vector meson differential cross 
sections is typically limited to \t\ < 1 GeV^ to avoid regions where large corrections were 
applied to account for the proton dissociation process. -B_d is related to the size of the inter- 
action region, as discussed below. It is sensitive to the gluonic proton shape and to the sizes 
of the dipoles contributing to the process. The dipole sizes depend, through the incoming 
and outgoing wave functions, on the value and the vector meson spatial extension. 

Fig. [6] shows Bd as a function of W for J/ip photoproduction. The average value of Bd 
from this process was used to fix the Bq parameter of the proton shape in the b-Sat model 
to Bg = 4.0 GeV~^. The model allows then to predict the expected B^s for other processes, 
like DVCS and and p production. The predictions are compared to data in Fig. [3 The 
value of Bd depends considerably on Q"^ and the type of process since the effective size of 
vector mesons and of the DVCS photon differs substantially. The dipole model reproduces 
the changes in the value of B^ very well. 

As seen in Fig. |6]the value of Bj:, in J/tp photoproduction shows some W dependence. 
This can be attributed to QCD evolution effects. In the DGLAP evolution scheme, due 
to strong ordering, the increase oi Br, with decreasing x ~ Mj, is expected to be 
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Figure 6: The t slope parameter Bd vs W for J/ip photoproduction. 



neghgible [38], as in the b-Sat model. In the BFKL evolution scheme the evolution in x 
could lead to a small increase of the size of the interaction area. This was taken into account 
in the b-CGC model, which is based on a phenomenological approach to the BFKL and BK 
equations [18-20]. Alternatively, the variation of with W can also be attributed to the 
increase of the contribution of the "pion cloud" at small x [39]. 

3 J/iIj as a probe of proton and nuclei 

The properties of diffractive processes described in the previous section single out J/ip pho- 
toproduction as an ideal probe for investigation of nuclear properties due to three main 
reasons: 

• The J/ip meson provides the smallest probe when compared to the other diffractive 
processes measured at HERA. This can be directly seen from the measurement of the 
sizes of the interaction region, B]j, shown in Fig. E] and [71 

• The observed number of well measured events is substantially larger for J/ip pho- 
toproduction than for other exclusive diffractive processes with similar probe sizej^ 

^ This is due to the fact that the cross section for the electroproducion process is approximately propor- 
tional to log((5^;j2,/Q^j„) ^'^d to the photoproduction cross section a'^P. In the ZEUS and HI measurements 
of J/4> photoproduction, ~ 10^^^ GeV^ and Q^ax ~ 1 GeV^. In this region the photoproduction 

cross section, a'^^, is almost constant. For the p, cj) or DVCS process one has to require that Qmin ~ 10 GeV^ 
to assure that the interaction is mediated by the small dipoles. In addition, the virtual photon-proton cross 
section is droping quickly with increasing Q^, a'' 'P 1/{M^ + Q^)", with n sa 3, which hmits effectively 
the range. 
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Figure 7: The t-slope parameter B^, for the DVCS process and for exclusive [35] and 
p [34, 36] meson production compared to predictions from the b-Sat and b-CGC models. 



In the central region of the HI and ZEUS detectors, the number of well measured 
7P ^ J/ipp ^ fifip and 7P ^ J/ip p ^ eep events is about 10 times higher than for 
7*p —^pp—^TTTTp and by a factor of 40 higher than for 7*p — 0p — KKp processes 
with + ]Vf2 > 10 GeV^. For the DVCS reaction, with > iq GeV^, this factor is 
around 100. 

• The J /if) meson decays with a probability of 12% into a leptonic pair, p^p^ or e+e^. 
These are very clean final states which emerge from quark annihilation. It can be well 
measured in the detector because the J/ip is a. very narrow resonance. The p meson 
decays mostly into a pair of pions which can be also well measured. However, this 
is a strong decay, the p width is large and the decay mechanism is more complicated 
than in the J/ip case. This may make this process less suitable as a probe of nuclear 
properties. 



3.1 Proton radius 

To show the potential of J/ip photoproduction as a probe of the properties of matter we 
discuss here a classic nuclear physics subject, the determination of the proton radius. 

The smallest sizes of the interaction region, Bn, are measured at HERA in the exclusive 
J/ip photoproduction process, in agreement with the dipole model expectations. The t 
distribution of this process is also one of the most precisely measured at HERA. In the 
following we will determine the proton radius using data of this process. 

The observed values oi B^ in the HI and ZEUS experiments show an increase with 
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growing energy, W. This increase is parametrized by the HERA experiments as 



Bd{W) = bo + 4a' • log(iy/90GeV) (13) 

For ZEUS 

60 = 4.15 ± 0.0512:?^ a' = 0.115 ± O.OlSlg;!]?^ 

and for HI 

60 = 4.63 ± O.Oel^i a' = 0.164 ± 0.028lJ]:[!^[] 

The W region extends from 30 to 170 GeV and the values of Bij increase with W by about 
25%, due to the QCD evolution effects. The proton radius should be related to the value 
of Bd without this effect, i.e. to the value of B^ at W near the threshold. However, it is 
not known whether the formulae flTSl) can be extrapolated below W = 30 GeV. Therefore 
we evaluate the proton radius from the lowest observed value of Bd at W = 30 GeV. This 
value of W corresponds also to the x value, x = 10~^, up to which the dipole picture was 
successfully tested. 

At = 30 GeV the size of the interaction region is Bd = 3.64 GeV^^ for ZEUS and 
Bo = 3.91 GeV^^ for HI. We can combine these values to 

Bn{W = 30 GeV) = 3.78 ± 0.3 GeV^^ 

As error we take the difference between the ZEUS and HI values. 

The proton size is related to the value of Bq which is smaller than Bd, due to the 
contribution of the size of the J/ijj. The difference at = 30 GeV is B^ — Bq = Ai? = 0.6 
GeV^^ [18] and is almost independent of W . The theoretical error can be estimated from 
evaluating Bq with two different wave functions as Ai? = ±0.2 GeV^^ [18]. 

The corresponding proton gaussian width is Bg = 3. 18 ±0.4 GeV~^, where we added the 
theoretical and experimental errors in quadrature. The transverse proton radius is then 



\/< 62 > = y y dWTGih) = ^1-Bg = 0.50 ± 0.03 fm 

The proton radius is usually determined from the electromagnetic charge form factor, 

GE{t) = l + ^<rl>t + Oit^) 

and 

\t=0 



^ dt 

This is the 3-dimensional proton radius, the transverse proton radius is given by 

dGE{t), 



dt 



|i=0 
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Therefore the transverse proton radius is related to the 3-dimensional proton radius by 
< bp >= 2 < Tp > /3 [39,40]. Guided by this we obtain the 3-dimensional proton radius as 
measured by the exclusive J/ip photoproduction 

^< r^g > = \/?>-Bg = 0.61 ± 0.04 fm 

We call this radius r2g to indicate that it is determined in the two gluon exchange process. 

The two gluon proton radius, r2g, is much smaller than the charged proton radius deter- 
mined from the electromagnetic form factor Ge, Tp = 0.875 ± 0.007 fm. One could argue 
that it is more appropriate to compare the two gluon proton radius to the proton radius 
which is determined from the Dirac form factor Fi, rp = 0.81 fm, instead of Ge, as proposed 
in ref [40]. This value is smaller than the standard charged proton radius bccaiise the Dirac 
form factor describes only the spin non-flip interactions (in the infinite momentum frame). 
Nevertheless, the spin preserving proton radius is still substantially larger than the two gluon 
proton radius. 

It is expected that the value of the proton radius is process dependent because the current 
which tests the proton has itself a structure which depends on its quantum numbers [41]. The 
smallest proton radius is determined from the axial form factor Ga measured in neutrino 
scattering up fip. It is called the axial radius and has the value of = 0.675 ± 0.02 
fm [41]. The smallness of this radius can be attributed to the fact that the axial current is 
not coupling to the pion cloud surrounding the bulk of the proton [39]. It is interesting to 
observe that the two gluon radius is still smaller than the axial radius in spite of the fact 
that the two gluon intermediate state couples as well to the bulk of the proton as to the 
surrounding pions. 



3.2 Nuclear structure 

The measurements of diffractive processes on nuclei could become an important source of 
information on nuclear structure and high density QCD. The interaction of a dipole with a 
nucleus can be viewed as a sum of dipole scatterings of the nucleons forming the nucleus. 
The size of the charmed dipole in elastic J/ip scattering is around 0.15 fm [17]. i.e., it is 
much smaller than the nucleon radius. It is therefore possible that dipolcs interact with 
smaller objects than nucleons; e.g., with constituent quarks or hot spots. Nevertheless, for 
the sake of example, we take here the conventional point of view and assume that the nucleus 
is build out of nucleons and that dipoles scatter on the ensemble of nucleons. We discuss 
the dipole- nucleus scattering with two examples to show the potential of possible nuclear 
investigations. To simplify the discussion multiple scattering effects will be ignored justified 
by the small size of the J/ip dipole, see also Section 3.3. 

In the first example we consider elastic scattering on the nucleus [17,21]. This is a coherent 
scattering process since the nucleus remains in its ground state. Ignoring possible multiple 
scattering effects, the dipole scattering amplitude for a given configuration of nucleons {bi}, 
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using eq. M and eq. [ID] is given by 



where A denotes the number of nucleons in the nucleus, Bp is the diffractive slope of the 
proton and ap is the total proton dipole cross section. The Fourier transform of this amplitude 
is 

i=l 

Here we changed the integration variable from h to h — hi while integrating each term of the 
sum. We obtain the matrix element for the elastic scattering (gg) + Aq ^ (qq) + Aq by 
averaging over all configurations of the nucleus ground stat^. 

- ^<^Ao = ^pe-^^-^'/2 Yl / d%..J'bAnoibi---bA)^Aoibi...bA) ■ e"^'''^ (16) 

i=i ^ 

We define the single nucleon distribution as 

j d%...d'bAd^^X{h...bA)'^Ao{bi...bA) = Ta(6i) (17) 

with normalization J d%iTA{bi) = 1. Because the wave function for protons and neutrons 
is completely antisymmetric and the difference between proton and neutron is presumably 
small we assume that Ta(6i) = T^(6i). The cross section for the coherent dipole scattering 
is then given by 

^ £j_^^-B,.A^ . d%TA{b)e-^'-'' , (18) 
dt IGvr J 

i.e. it is proportional to the square of the Fourier transform of the single nucleon distribution. 

The cross sections for coherent J/ip scattering on nuclei, computed in the dipole model 
by averaging the dipole cross section, eq. [HI over the incoming photon and outgoing J/ip 
wave functions, are shown in Fig. [8] as a function of t. Note that in the nuclear case, the size 
of the vector meson can be ignored. We assumed here that the single nucleon distribution 
can be identified with the Woods-Saxon distribution, see Appendix. The diffractive slope 
at t = depends on the size of the system. The figure shows for small t ~ 1/-Ra ^ very 
steep t dependence, ~ exp(— ti?^/3), and then several diffractive minima. For larger nuclei 
the Woods-Saxon shapes are approximately similar to a box with a size Ra- The Fourier 
transform of a box is given by the Bessel function J\{Ra ■ A) which has zeros at Ra ■ A = 
3.8, 7.0, 10.2.. .M This gives the approximate positions of the diffractive minima for the 
exclusive J/ip photoproduction on Calcium (A=40) and Gold (A=197) seen in Fig. [Hi The 
parameters of the Woods-Saxon distribution, nuclear radius and skin depth, were determined 



^To simplify the notation wc write ^'(■•■) as a function of the transverse variables, bi, only; i.e. we assume 
that the longitudinal dimensions, Zi, are already integrated out 
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This is like the Frauenhofer diffraction on a circular aperture. 
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mainly by scattering of the charged matter and can be pretty different when measured in 
dipole interactions. We therefore plotted in Fig. M the predictions computed with slightly 
altered values of these parameters. 

In the second example we evaluate incoherent dipole scattering [46], [13,21,22]. We start 
with the quasi-elastic scattering of a dipole on a nucleus 



(qq) Ao ^{qq) A„, 

n 

where An can be either a ground state or any excited nuclear state or any breakup of the 
nucleus into nucleons or nucleonic fragments. Pion and other hadronic production is not 
allowed. Like in the elastic case, the matrix element for transition to a state An is 

^ r 

- = cxpe-^-^'/' E / d%..AnAblSA)^A,{bl..bA) ■ e-*-^ (19) 

i=l 

The quasi-elastic dipole cross section 

dt 167r ^ I I 

n 

o A A 



3^ 

167r 

« 3 



n 

can be evaluated using the completeness relation 

^aM---VaWaSKSa) = 6{h - b\)...6{bA - 9a) (20) 



giving 

E^^^ = l^e"''^^'EE / d%,..A%AnSbi---bA)-^AAhi.SA)-e-^^^^^^^^ (21) 

71 i j 

Defining the two nucleon distribution as 

J d%,..A%Ad^^\^{K-iA)-^AAbi---bA) = Tf (&i,&2). (22) 
with normalization j d%id%2TA^ i^i,b2) = 1, we obtain 



^ — = ^e-^-^' A^A{A-\) / d%d%T^^\b,, h) ■ e-(^^-^^)-^ 



(23) 



The first term in the square brackets, proportional to A, emerges from the summation of 
terms with i = j in eq. EH The second term, proportional to A{A — 1), is obtained using the 
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antisymmetry property of the wave function and the assumption that the difference between 
protons and neutrons can be neglected, T^\bi, 62) = T^\hi, bj) for i ^ j. 

The incoherent dipole cross section is obtained by subtracting the ground state contri- 
bution from the sum of eq. 



^^S^ = ^e-^^"-' I d%d% {a (t^(6i)T^(62) -Tf (&i,&2)e-^(^^-^^^)-^) 

+ A' (Tf (61,62) -T^(&i)T^(&2)) e-'^'^-''>^] (24) 



dt 16n 



When the momentum transfer A ^ the incoherent dipole cross section, eq. [211 goes 
to zero reflecting the fact that without a transfer of transverse momentum to the nu- 
cleus excited states cannot be produced. At larger values of the momentum transfer, 
|A| > 200 MeV, the contribution of the second term in the upper bracket of eq. [2U 

/ d%d%A (^Tj^^(6i,62)e-*(^i-''2)-^j, starts to be relatively smaller than that of the first 

term because the oscillatory factor suppresses the contributions with |6i — 62I > 1 fm. The 
contribution of the second term grows then as A^^^ and can be neglected for larger nuclei 
compared to the contribution of the first term which grows as A. Thus when |A| > 200 MeV 



E 



dt dt 

+ ^e-^^^^A^ ^^61^^62 (rf (61, 62) - r^(6i)r^(62)) e-^^'^~'^>^ (25) 

and the deviation from the elastic nucleon cross section result is due to two-body correlations 
only. 

Experimentally we expect to be able to distinguish cases where the nucleus remains intact 
and cases where the nucleus breaks up. In the nuclear breakup process, there are around 
0.3v^ free neutrons and 0.2v^ protons in the final state [47], as well as other fragments. 
These particles and fragments have high momenta and different charge to mass ratios than 
the nuclear beam, and should therefore be measurable in specialized detectors. However, we 
do not have a one-to-one correspondence between an intact nucleus and a coherent scattering 
process since incoherent processes can lead to an intact nucleus in an excited state. It 
remains to be determined how well excited states of the nucleus can be identified, measured 
and possibly statistically subtracted to extract the fully separated coherent and incoherent 
processes. 



3.3 Saturation 

One of the main results of HERA is the observation that the gluon density increases quickly 
with decreasing x. This suggests that at some small x the dipole could undergo multiple 
interactions. The degree of saturation is characterized by the size of the dipole rs which, at 
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Figure 8: The prediction of the dipole model for the t distribution of coherent J/ip pho- 
toproduction on nuclei assuming that the single nucleon distribution can be identified with 
the Woods-Saxon distribution given in Appendix. The statistical errors of the simulated 
measurements are based on the assumed collected sample of 10^ events. The upper dashed 
line shows the sum of the coherent and incoherent process in case of no correlations. 
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Figure 9: The distribution of the J/ip dipole sizes; Left: as given by the 'jJ/4' overlap 
function; Right: as given by the photoproduction amphtude. 



a given x, starts to interact multiple times. The dipole size rs is defined, by convention [18], 
via the relatior0 

daaaix.rQ.h) , ^ 

=21- exp -1/2 ^ 0.8. 26 

The saturation scale is defined then as Q| = 2/r| and is a function of x. A high value of the 
saturation scale means that gluon density is so high that even small dipoles interact many 
times. 

In various analyses of HERA data the saturation scale, in the proton center, was found 
to be Ql ~ 0.5 GeV^ at x ~ 10~^, i.e. in the EIC range. In the LHeC range, which extends 
to a; ~ 10~^, the saturation scale could reach Ql ^ 2 GeV^ [18]. The saturation scale deter- 
mined in the inclusive 7*p reaction should be compared to the scale of J/ijj photoproduction 
given by the effective size of the meson in this reaction. This size is determined by the over- 
lap of the photon and J/ ip wave functions and by the dynamical effects. Fig. ^^eft shows the 
distribution of the dipole sizes selected by the overlap of the photon and J/ ip wave functions 
defined as r J The median value of this distribution is r^^^ = 0.7 GeV~^, i.e. 

around 0.15 fm. The contribution of very small dipoles, in the jp reaction, is suppressed by 
the dynamic of the reaction, because the dipole cross section is proportional to at small 
r, eq. [91 The distribution of the dipole sizes selected by the amplitude, eqJTl is shown in 
Fig.^^ight. The median value of this distribution, r^g^, was found to be between 1.2 and 
1.4 GeV~^, depending somewhat on x. In the jp reaction the effective scale characteristic 
for J/ip photoproduction is then Qlff = '^/ijmed)'^ ~ 1 — 1.5 GeV^, depending on x. The 

From the unitarity limit the highest value of a dipole cross section \s da / (Ph = 2, see eq. [TTJ 
^^The evaluation was performed in the b-Sat model with the' Boosted Gaussian' wave function , at i = 
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Figure 10: The survival probability S"^ as a function of x for various dipole sizes. 

saturation scale at ETC is sizably smaller then Qljf, so that saturation effects are not ex- 
pected to be large. However, the measurement on nuclei could enhance the saturation scale 
substantially, as discussed in ref. [21]. Therefore, the high precision which can be achieved 
in the measurement of J/ip exclusive photoproduction makes this process interesting as an 
alternative testbed of saturation effects at EIC. At LHeC saturation effects should be clearly 
visible in the scattering on proton and presumably very strong in nuclear reactions. 

In ref. [9] it was proposed to investigate the effects of saturation in a systematic way by 
extracting from data 5*^(6), the square of the 5* matrix. The 5* matrix is directly connected to 
the dipole cross section by dagg/cPb = 2[1 — ^S{b)] as shown in eq. [51 The S matrix can also 
be used to define the saturation condition because its square has a meaning of the survival 
probability; i.e., the probability for no interaction. The equivalent definition to eq. [261 then 
reads S"^ = e^^ ~ 0.37. Fig. [10] shows the survival probability for dipoles of different sizes 
scattering on the proton, at 6 = 0, where the gluon density reaches its maximum [17]. 

To reconstruct the S matrix at 6 = for protons it is necessary to measure the t distri- 
bution up to about |t| 2 GeV^ 0. At HERA this measurement could not be performed 
because the measurement of the t distribution had large systematic errors for |t| > 1 GeV^ 
and a low statistical significance. The systematic uncertainties which plagued HERA ex- 
periments can be avoided in the specially designed experiment discussed below. In addition 
it should be possible to increase the collected number of events by around two orders of 
magnitude in comparison to HERA which should allow to reach the |t| = 2 GeV^ region in 
the proton measurement. 



In case of nuclei the measurement will require a much smaller t range than in the proton case. 
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4 Production Cross Section 



In this section, we review the production cross section for the photoproduction of J/ip and 
discuss the kinematic range where high precision measurements should be possible. The 
ep — > epJ/ip cross section can be written as 



a 



1 



where 
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and 1/ is the inelasticity. We are interested in very small Q^, e.g., < 10~^ GeV^, and we 
can assume that in this range , the photoproduction cross section is independent of Q^. We 
also assume that the longitudinal cross section is negligible since we are dealing with almost 
real photons. Then integration gives 
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Prom the ZEUS data [33], we have 
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Writing this in terms of y, we have 
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where W is in GeV and s is in GeV^. This form is valid at the large W measured at HERA. 
An extrapolation to lower values of W is in agreement with measurements performed by 
E401 [50], so we will assume this form for all W. 

We consider two different values of Q'^ax- 1^ case, we assume that we do not measure 
the scattered electron precisely, and we can only limit the Q^. For this case, we take Q'^^x — 
10"*^ GeV^ so that the maximum from the electron is 10 McV. In the second case, we 
assume we do measure the scattered electron well, so we can afford to go to higher Q^, and 
we take Ql^ax = 10~^ GeV^. In the Q'^ range considered here, the scattered electron energy 
is given by E' ^ {I - y)E^. 

High precision measurements of the J/-^ decay products are assumed to be made in a 
central detector, such as a thin-walled time projection chamber, located in a strong solenoidal 
magnetic field. The detector parameters are discussed in the next section. To measure the 
J/ip decay products with good acceptance in the central detector, we require that the J/ip 
has limited boost. For simplicity, we take = t = to see the limits on y [49]. Conservation 
of energy gives 
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where E^, Ep are the incoming beam energies, and E'^, Ey are the scattered electron and 
vector meson energies. This can be rewritten as 

Ev - yEe 



X 



E„ 



We also have the requirement: 



My — {xp + qf ~ 2xp ■ q ~ sxy 



or 



X 



Ml 
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Putting these expressions together gives 
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where My is the mass of the J/V'- This leads to the following constraints: 
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We now integrate the differential cross section given above in this W'^ range, with the 
result: 



cr (visible) = A 
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where we have 



A = 75. nb 
b - 0.35 

Y ^max 
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To see the effect of the ccntrality requirement, we require that the momentum of the 
J/ if} is less than 4 GeV, so that the decay particles are not too boosted. We performed a 
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scan over electron beam energy, keeping the proton beam energy fixed at Ep = 100 GeV. 
The accepted y ranges are shown in Figure [TH The results for the visible cross section (not 
including branching ratio) are given in Fig. [121 Adding the branching ratio of about 6 % 
implies that 10 fb^^ would give 10^ well reconstructed J/tp —>■ ii^ ii~ events in the central 
part of the detector. 




Figure 11: The y range with good acceptance for reconstructing J/ip decays in the central 
region is plotted as a function of the electron beam energy (shaded area) at Ep = 100 GeV. 
The two curves indicate the minimum y values for detecting the beam electron scattered at 
degrees assuming (upper curve) that a 1 Tm dipole is placed near the IP and that the 
scattered electron should be at least 5 cm from the beam 5 m downstream of the dipole, 
and (lower curve) that a 5 Tm dipole is placed near the IP and that the scattered electron 
should be at least 5 cm from the beam 10 m downstream of the dipole. 



5 Experimental Considerations 

In the following, we consider the requirements on the beam momentum spread and the 
detector requirements for the scattered electron, the measurement of the J/ip, and the mea- 
surement of the scattered nucleus or proton. As discussed in Sections 2.2 and 2.3 the proton 
or nuclear shapes are measured by the deflection of the J/ ip from its original direction. The 
t measurement will rely on an accurate determination of the transverse momentum given to 
the J/ip since 

t=iP^- Pff ~ - {{^P.f + {^Pyf) 

where PiiPf are the momenta of the incoming photon and outgoing J/ip, and Apx,Apy are 
the changes in the transverse momentum components during the scattering process. As is 
clear from these expressions, we require knowledge of 

1. the outgoing J/tp momentum, and 
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Figure 12: Visible cross section as a function of the electron beam energy at Ep — 100 
GeV. 

2. the incoming and outgoing electron momenta since 

Pe,i - Pe,f - Pj/^ ^Pi-Pf ■ 

Usually, t is viewed as a change of the momentum of the incoming and the outgoing ions 

or protons, t = —{p7.i — pTjY ■ The measurement of the ion momenta is technically very 
difficult. Here, it is not necessary because p/ j ^ Pi'j = Pe,f '^Pj/ij} ~ Pe,i- We, therefore, do 
not foresee a precision measurement of the ion p^- 

Measuring the outgoing J/ip and electron (or guaranteeing that the outgoing electron 
has very small transverse momentum) then allows a measurement of t. Nevertheless, the 
forward scattering region for the ions will require substantial instrumentation in order to 

1. guarantee that we have an elastic scattering event. 

2. measure the correlation between the variable t and the transverse momenta of debris 
from the nuclear dissociation in inelastic events. 

We discuss the separate parts of the detector system according to their function. Since we 
are deahng with photoproduction, the large majority of the events are concentrated close to 

= GeV^. The scattered electron will therefore have a median transverse momentum of 
only a few MeV. As discussed below, it should be possible to measure this scattered electron 
over the bulk of the interesting kinematic range. Even in cases where the electron is not 
detected, it will typically only carry a very small transverse momentum. The measurement 
of the t distribution will therefore rely almost exclusively on a precise measurement of the 
J/ip decay products. 
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5.1 Scattered electron 



The scattered electron will be at a very small angle to the direction of the electron beam, but 
can be pulled out of the beam if it has lost enough energy. There will be a minimum value 
of y where it is possible to see the scattered electron even if the electron is scattered at 0°. 
For this purpose, we imagine that there will be a large-aperture dipole magnet placed near 
the beam interaction point which will bend the electrons (and can also act as a separator 
of the ion and electron beams). The transverse separation between a beam electron and a 
forward scattered electron can be written as: 



A = D 

where 



' . . cos ^0 . n 
sm Uy sm 

cos 6y 



^ 0.3-B-L 

sm = 

P 

with the dipole field strength B measured in T, the length of the magnet L in meters and the 
momentum of the electron p in GeV/c. Oq is used to denote the scattering angle of the beam 
electrons after passage through the dipole, and 6y is the scattering angle for momentum 
p = {1 — y)Ee. The distance D is measured from the center of the dipole. If detectors 
are placed within a few cm from the electron beam a few meters downstream of the dipole 
magnet, then scattered electrons at not too low y can be measured. The minimum value 
of y required for acceptance of such a system depends on the strength of the dipole field, 
the distance of approach to the beamline, and the drift distance D. Figure [TT] indicates the 
minimum y values as a function of the electron beam energy for two representative sets of 
parameter values. We conclude that measuring the scattered electron will be feasible over 
the bulk of the interesting y range provided the interaction region allows for such a dipole 
magnet and detector arrangement. 



5.2 J/i/j reconstruction 

We focus on the measurement of the J/ip via the decay into muons, J/ifj ^ . We 

assume no particle identification will be necessary, since it will be easy to identify the J/ip 
via the invariant mass of the reconstructed state. A continuous background from the Bethe- 
Heitler process eA — > eAfi~^fi~ will have to be handled in the analysis of the data. The 
expected resolution of the drift chamber can be estimated from the measurement precision 
term 

(a /m) - ^'""'^ 

\^Pt I ft J meas 



0.3L^B V iV + 4 
and the multiple scttering contribution 



KM)m5 = |||^1.43A[i + 0.038 log(L/Xo)] 

as 

(^Pt/Pt = {(yptlPt)meas © (f^ptM)AfS- 
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Here B is the magnetic field in Tesla, L is the lever arm in meters, ar4, is the spatial resolution 
in meters for a single point, Xq is the radiation length in meter, N is the number of points 
and f3 is the velocity of the particle. 

In the J/ip photoproduction process at small t the transverse momenta of the decay 
muons range from very small values up to a maximum of about 2 GeV/c. In this momentum 
range, multiple scattering is a critical issue for momentum resolution. We therefore envision 
a TP C- type detector with a thin inner wall as the central tracking detector. Assuming the 
following parameters 

1. outer radius R = 2 m 

2. solenoidal field B ^ 3.5 T 

3. gas density Xq — 450 m 

4. point resolution cr = 100 /xm 

5. measurement N — 200 points. 

yields a track momentum resolution [48] 

apjpt = 0.005 • pt e 0.045//3 % 

which will give a pt resolution for the J/ ip of typically < 1 MeV. 

The muons would not escape such a detector in the radial direction. An electromagnetic 
calorimeter placed outside the TPC could be used to reject or measure radiated photons from 
the interaction vertex and from passage of the muons through the detector. The detector 
design could naturally be extended to include hadronic calorimetry, particularly in the ion 
direction, resulting in a general purpose detector capable of high precision measurements for 
all types of exclusive processes as well as inclusive cross section measurements. 

5.3 Ion forward direction 

The main requirement for the instrumentation in the forward direction of the ion beam is 
that dissociative and inelastic events be rejected with ~ 100 % efficiency. We envisage again 
a dipole magnet which can be used to separate neutral particles as well as nuclear fragments 
from the main beam. This dipole could also serve as a beam separator to guide the electrons 
and ions to their individual beampipes. Ion dissociation will produce different types of 
fragments, including neutrons and charged ions with different charge/mass ratios than the 
beam ions. Neutrons can be recorded in a calorimeter at zero degrees to the beamline at the 
interaction point, and located many meters downstream of the interaction point. Fragments 
with different charge/mass ratios will have different deffection angles in the dipole than 
the main beam, and can therefore be measured with tracking detectors placed close to the 
beamline. In inelastic scattering events, particles with opposite charge to the beam ions can 
be produced, and these interactions can be vetoed using tracking detectors or calorimeters 
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placed behind the dipole. We therefore expect that for large nuclei, rejecting non-elastic 
events will be rather straightforward since typically several neutrons and charged particles 
will be produced. 



5.4 Beam requirements 

The experiment discussed requires substantial instrumentation extending to many meters on 
either side of the interaction point. It will be critical to avoid 'dead zones', where scattered 
particles could escape detection. This will place severe restrictions on the accelerator design, 
and will naturally lead to limitations in luminosity. The detector design foresees a central 
region with a strong solenoidal field, and dipole magnets on either side of the interaction point 
(perhaps +-2 meters away). The design of these magnets will necessarily have to be done in 
conjunction with the accelerator group. Issues such as synchroton radiation loads will need 
to be evaluated. In addition, it is critical that the electron beam transverse momentum be 
limited to a few MeV, since this quantity cannot be measured but enters in the calculations 
of t. 



6 Summary and Discussion 

We discussed here the physics potential of the t measurements. At HERA, the measure- 
ment of the t distributions allows to determine the two gluon proton radius. This radius is 
substantially smaller than the proton radius determined in the electromagnetic interactions. 

To exemplify the potential of such measurements on the future electron-ion colliders we 
considered two simple examples of coherent and incoherent J/ip dipole scattering on nuclei. 
In the coherent case the nucleus remains intact while in the incoherent case the nucleus 
goes into any excited state or it disintegrates into nucleonic fragments or nucleons without 
production of additional hadrons. 

In the coherent case we computed the predictions of the dipole model assuming that nu- 
cleons are distributed within the nucleus according to the Woods-Saxon distribution. Fig. [HI 
which is mainly determined by scattering on the electric charges. It is well possible that, as 
in the proton case, the nuclear radius and skin parameters can be quite different when mea- 
sured by elastic J/ip scattering. The shape of the expected cross section and the differences 
between different predictions indicate that a measurement resolution of around 10 MeV for 
the pt of J/ip should be sufficient. This number emerges from the requirement that the first 
diffractive minimum should be properly resolved. 

Incoherent J / ip scattering is equally interesting since the very good pt resolution com- 
bined with the full acceptance detector discussed here should allow a systematic study of 
the two body correlations. The measurement of the long range nuclear correlations could 
require a measurement resolution of 0(1) MeV in px of the J/ip decay products. 
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The t distribution of the nuclear breakup process will be measured in correlation with 
the number and momenta of the breakup protons and neutrons. This should allow to study 
the dissociation process as a function of the transverse momentum transfered to the nucleus. 

Although we concentrated here on the measurement of J/ip decays it should be clear that 
the detector which is optimal for measurements of J/ ip elastic scattering will also be optimal 
for measurements of other interactions, like light vector meson production and inclusive 
diffractive or inclusive total cross sections. The quality of the inclusive measurements will also 
profit from the large acceptance detector discussed here because the coverage of the almost 
entire rapidity range for charge and neutral particles will reduce the systematic uncertainties 
in the F2 and F2 measurements. 

The nuclear effects will also be seen through the measurement of the absolute value of 
cross sections. The dipole model predicts these values precisely provided the scattering takes 
place on nucleons within the nucleus which have the same properties as free protons. Any 
deviation from the expected value carries information about nuclear effects. For example, 
the total dipole proton cross section, dp, can be different for a nucleon in a nucleus than 
for a proton because a nucleon within a nucleus can have a different size than a free proton 
or neutron. This would change sizably the value of the nuclear dipole cross section and 
therefore also the values of the observed diffractive cross sections. By the same argument 
the measurement of F2 on nuclei is also determined by the nuclear properties and will lead 
to very interesting saturation effects as discussed in [17,21]. Saturation is dependent on the 
size of the scattering objects, therefore its measurement and the absolute value of the cross 
section could indicate on what objects the scattering takes place. 

The measurement of the transverse shape and correlations combined with the measure- 
ments of the total cross sections and the shadowing effects should allow to determine the 
inner structure of gluonic fields which keep the nuclei together. 



7 Conclusion 

Wc have shown that dipole interactions of small dipoles with nuclei can become an impor- 
tant source of information on the nuclear structure. At high energies, the dipoles interact 
with nuclei by a well understood QCD process in which two gluons, with large transverse 
momenta, are exchanged. The difference of the gluon momenta can be precisely determined 
by measuring the transverse momentum of the elastically scattered vector meson. The J/%1) 
photo-production process is particularly well suited to perform this measurement because 
it has the largest cross section and the best measurement precision. The Pt of J/"^ can be 
measured with high efficiency and a precision of 0(1) MeV using presently available tech- 
niques. This allows to determine the t distributions with a precision 0(1) MeV^ starting at 
i fa 0. The upper range of the available t values depends on the process and can reach 0(3) 
GeV^ with the luminosity foreseen at ETC or LHcC. This allows to investigate large and/or 
small gluonic structures which keeps the matter together. It is worth emphasizing that large 
structures are investigated by a highly virtual, well understood interaction because a small 
Pt of the meson is a result of the difference of two large transverse gluon momenta. 
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Therefore, dipole measurements will provide high quality data of basic nature which could 
lead to solution of the long standing puzzle; how strong interactions form matter. 
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9 Appendix 



9.1 Derivation of the dipole representation 




Figure 13: Elastic photon-proton scattering at t = 0. 



We give here a derivation of the dipole representation for the elastic 7*p — > 7*p scattering 
to exemplify the main properties of the dipole picture. Dipoles incorporate naturally the 
interference effects between different quark-gluon couplings which are essential to obtain a 
proper description of DIS reactions at low x, namely colour transparency!^ . For the sake 
of example let us consider the scattering of a longitudinally polarized virtual photon on a 
proton. Following the derivation given in [14] and [8, 15] we write the total cross section in 
the kT factorization form according to the Feynman diagrams shown in Fig. 1 as 

- ^"-^Y.^)jfM)j^ij<^^^Q^^^l - (4 - ^) (27, 

"'^^ Although we want to limit ourselves in this paper to ep and eA scattering only, let us note that the 
dipole concept, because of its importance in understanding the QCD evolution [23] and its success at HERA, 
is now extensively investigated also in the context of pp scattering at the LHC [24] . 
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where D{k) = k"^ + Q"^ and = z{l — z)Q^ + rrij. nif denotes here the flavour dependent 
quark mass and f{x, t^) denotes the unintegrated gluon density. In the 7*p colhnear frame k 
and / are the two-dimensional transverse momentum vectors of the quark and the exchanged 
gluon and z and {1 — z) are the fractions of the light-cone momentum of the photon carried 
by the quarks. 

In the elastic forward scattering there are just two possibilities for gluons to couple to 
the quarks shown by the two diagrams of Figure [131 Both diagrams should be taken into 
account because only then the important property of colour transparency emerges from the 
cancellation between the two propagators in eq. [271 In the limit k ^ / the cancellation leads 
to (Tl ~ 0, reflecting the intuitively clear fact that a gluon cannot see a quark pair when its 
wave length is much larger than the distance between the quarks. In this case, the quark 
pair will appear as neutral to the gluon. 

Colour transparency emerges more naturally when the transverse quark momentum k is 
replaced by its Fourier conjugate r, the transverse separation between the two quarks. We 
have then 

d?k ^1 
— exp(ifcf)-— ^ = Ko{r ■ Q) 
27r D{k) 



and 



exp(— zA;r)(l — exp(— z/r))/ro(^ ' Q) 



D{k) D{k + l)J J 27r 
with Kq being the Bessel-Mc Donald function. 

The change of variables from k to f leads to the dipole representation, 

aL = ImA-^'Pix, g, A = 0) = J d^r J dzj^ Q>gg{x, r)^i(r, z, Q^) (28) 

in which the wave function ^> l describe the probability amplitude to find a qq pair within a 
virtual incoming or outgoing photon. 

^l=^-^e)AQ^z\l-zfKl{r.Q) 

The dipole cross section cjgg describe the interaction of the qq pair with the proton mediated 
by the two gluon exchange. 

= f / - e-'0(l - = ^ / - Jo(/r)) 

By introducing the relation between the integrated and unintegrated gluon density 

xg{x,Q^)= r dff{x,f)/f, (29) 
Jo 

and assuming that as depends on the dipole size only, the dipole cross section can be further 
simplified by approximating (1 — Jo{lr)) ^ (/r)^/4, which is valid for P < 

aqq = —as{l/r^)r^xg{x, 1/r^). (30) 
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Here xg{x, /i^) is the gluon distribution which evolves in /x^ = according to the DGLAP 
evolution equation. 

In this representation colour transparency becomes a property of the dipole cross section, 
for small dipoles r ^ also agq — > 0, as is intuitively clear. The last form of the dipole cross 
section was first derived in ref [16], in an alternative way. 



9.2 Woods-Saxon Distribution 



The distribution of nucleons in the nucleus pA{r) is usually parametrized according to the 
Woods-Saxon distribution [42] 

Pwsir) = . . , (31) 

exp (^^^ j + 1 

with 6 = 0.54 fm, Ra = (1.12 fm) A^^^ — (0.86 fm) A"^/^ and N is adjusted to normalize the 
distribution to one 

' d'fpws{r) = l. (32) 



The transverse distribution is defined as 



/ + C!0 
dzpwsiVb^T^). (33) 
-co 

Figure [14] compares the proton shape Tp{b) with the transverse density Tws{b) for several 
light nuclei. 
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